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KINETIC  THEORY  FOR  ELECTROSTATIC  WAVES  DUE  TO 
TRANSVERSE  VELOCITY  SHEARS 


INTRODUCTION 

..  \ 

-^hear  in  the  flow  velocity  of  a  fluid  leads  to  the  low  frequency  and 
long  wavelength  Kelvin-Helmholtz  (K-H)  instability  T"~-  The  velocity  shear 
can  be  generated  in  a  number  of  ways.  In  a  plasma  the  existence  of  an 
inhomogeneous  electric  field  component  transverse  to  the  ambient  uniform 
magnetic  field  can  provide  a  transverse  velocity  shear.  The  evolution  of 
the  K-H  instability  in  this  configuration  has  been  extensively  studied^'1^ 
Recently  some  space  observations  '  and  laboratory  experiments  seem  to 
indicate  that  ion-cyclotron-like  waves  are  observed  for  subcritical  field 
aligned  currents  and  therefore  the  origin  of  these  waves  are  somewhat 
mysterious.  A  crucial  feature  of  these  observations  and  experiments  was 
the  presence  of  a  transverse  component  of  a  zero  order  electric  field.  In 
order  to  study  the  role  of  the  transverse  electric  fields  in  the  generation 
of  the  ion-cyclotron-like  waveSj ’  we  suggested  a  mechanism  based  on  the 
coupling  of  the  negative  energy  ion  Bernstein  modes  (or  the  ion  cyclotron 
modes)  in  the  region  where  the  d.c.  electric  field  is  localized,  with  the 
positive  energy  ion  Bernstein  modes  (or  the  cyclotron  modes)  in  the  region 

S5— - 


where  the  d.c.  electric  field  is  absent 
energy  wave  growth  in  an  inhomogeneous 


This  is  similar  to  the  negative 
6 


mirror  geometry 


In  our  initial 


theory  we  idealized  a  typical  electric^  field  profile  by  a  piecewise 
continuous  function  for  simplicity.  Tire  gradients  of  the  electric  field 
were  ignored  so  as  to  avoid  the  K-H  /(nodes  for  which  the  second  derivative 

,  0, 


of  the  electric  field  is  necessary. 


’Here  we  use  kinetic  theory  to  obtain 


the  general  dispersion  relation  rigorously,  for  the  electrostatic 
oscillations  in  a  plasma,  in  the  form  of  an  integral  equation  for  an 
arbitrary  electric  profile.^  In  various  limits  we  reduce  the  integral 


\ 


equation  to  second  order  differential  equations  to  obtain  the  eigenvalues. 

The  integral  equation  will  be  solved  in  a  subsequent  paper. 

Manuscript  approved  October  2,  1987. 


THEORY 


The  equation  of  motion  of  a  charged  particle  in  a  uniform  magnetic 
field  in  the  z  direction  and  a  nonuniform  electric  field  in  the  x  direction 
is  given  by 


d-E  =  £  E(x)x  +  a  V  X  z,  (1) 

dt 

where  2  =  eBQ/mc  is  the  gyro frequency,  e,  m  and  Bq  are  the  charge,  mass  and 
the  ambient  uniform  magnetic  field.  The  constants  of  the  motion  are  (i)  H 
=  (v2  +  v2+  vz^)/2  +  eY(x)/m,  the  total  energy,  where  E(x)  =  -  a^xj/ax; 
and  (ii),  =  x  +  v^/S,  which  is  obtained  by  integrating  the  y  component 

of  (1).  Using  v  =  2(X  -  x)  in  the  expression  for  H  we  obtain  the 

y  & 

Hamiltonian  for  an  equivalent  one  dimensional  problem, 


H  =  v2/2  +  22«  -  x)2/2  +  J  Y  (x)  .  (2) 

X  5  nl 

Minimizing  the  potential  of  (2)  we  obtain  the  guiding  center  position 


S  =  x  + 


v 


X 


-  V*> 

2 


(3) 


another  constant  of  the  motion  which  is  an  implicit  function  of  X  and 

g 

therefore  is  not  an  independent  constant  of  motion.  Here  VD( £)=-cE( £)/B  . 

In  order  to  recover  the  fluid  K-H  modes  as  the  fluid  limit  of  the 
kinetic  formalism,  we  will  need  to  construct  an  equilibrium  distribution 
function  using  the  constants  of  motion,  such  that  the  equilibrium  density 
is  uniform  or  nearly  so.  However,  we  would  also  like  to  be  able  to  study 
the  more  general  case  of  an  equilibrium  with  an  arbitrary  density  profile. 
Therefore,  we  will  choose  the  distribution  function  to  be  of  the  form 


constant.  Then  we  obtain  an 


f0U,H)  =  nna)Fna,H),  such  that  /FndJv  = 


0 


0" 


appropriate  fg  for  the  study  of  the  classical  K-H  mode  by  setting  ng(£)  = 
constant,  and  an  appropriate  f^  for  the  general  case  by  relaxing  this 
condition.  Such  a  distribution  function  which  leads  to  an  equilibrium 
density  uniform  to  0(c)  for  a  constant  n^,  where  c  (=p^/L)  is  the  smallness 
parameter,  pi  (=vt^/S^)  *s  t*le  ion  gyroradius  and  L  is  the  characteristic 
length  associated  with  the  external  electric  field  can  be  found  by  a 
systematic  procedure  and  is  given  by, 


ioa,  H  )  =  N  exp  [-0H  ]g(0, 


(4) 


3/2  2 

where  N  =  no(g/2ix)  ,  g  =  1/v^.,  v^  is  the  thermal  velocity  and, 


g(0  =  exp  [e{eYU)/m  +  V*(0/2}]n<0 


-1/2 


where  h(  O  =  1  +  VE(f;)/2- 


The  quantity  t]  parameterizes  the  magnitude  of 
the  velocity  shear.  Note  that  there  are  two  crucial  parameters  in  this 
problem;  (i)  t|  and  (ii)  e.  We  will  allow  h  to  be  arbitrary  but  positive 
while  assuming  e  <<  1.  The  equilibrium  distribution  can  be  expressed  as, 


N  exp(-  |  w^)exp(-  |  v^) 


(5) 


/h(0 

where  we  have  expanded  the  x  dependence  of  (4)  around  £,  and  neglected  terms 

3  2 

of  0(e  )  and  higher.  Here  is, 


9  9  9  ^  t-i  (  )  (  O  U  \ 

W.  =  V  +h((;)u2  +  — ^ -  fu  <U  >  -  )  . 

1  X  IVS/  y  _2  l  y  y  3  ) 


(6) 


where  uy  =  vy  “  ^vy^  an<*  indicates  time  average. 

Integrating  (5)  over  all  velocities  we  can  show  that  the  equilibrium 

2 

density  distribution  n  =  Oq{  1  +  o(e  )  },  is  uniform  to  order  e.  It  is 

possible  to  devise  a  distribution  function  with  density  uniform  to  any 
desired  higher  order  in  e,  but  this  is  not  necessary  here.  For  generality, 
in  the  following  we  shall  consider  a  nonuniform  equilibrium  density 
profile,  i.e.  n^  =  nQ<E,). 

Now  using  the  definitions, 


<j>(r',  t')  =  exp  {-i(«t'  -  kyy')}<|>(x') 


<Kx')  =  dk/  exp{ik^x'  )<Kkp  , 


where  <|>  is  the  electrostatic  potential  for  the  perturbed  electric  field  and 
linearizing  the  Vlasov  equation,  we  obtain  the  perturbed  distribution 
function 


fx(x,  v)  = 


0mfo[Idkx  exP  f 1  <k^x) }  <J>(k^)  +  ijdk^<J>(k^)(w  -  kyVg)Jdt'A  j  ,  (8) 


where,  x  =  t '  -  t, 


A(t',x',y',z')  =  exp  ji(k^  x'  +  ky(y'  -  y)  +  kz(z'-z)  -  wx) 


1  Un  ^  pQ 

vga)  =  enuTTo  aT  =  ve(U  "  2vi  +  — ~ 
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Here  cn=p/Ln,  and  Ln  =  ((dnQ(S)/d  0/^(0)  *  is  the  scale  length 
associated  with  the  equilibrium  density  profile.  In  order  to  evaluate  the 
time  integral  in  (8)  we  need  the  orbits.  The  x  component  of  (1)  gives  the 
net  force  in  this  direction.  Since  v  is  oscillatory  in  equilibrium,  <v  > 
=  0.  This  leads  to  the  expression  for  an  average  equilibrium  drift  <v^>  = 
<V„(x)>  in  the  y  direction.  Now  expanding  V_  around  £  we  obtain, 

h  fcj 


<v  > 

y 


VF(S)  +  -  <  (v  -  V„(0)  >  ♦  0 (eJ)  . 

E  2Q2h(0  y  E 


(11) 


Transforming  the  equation  of  motion  (1)  into  a  frame  moving  with  <v^>  1°  y 
direction  (i.e.,  Vy  Uy  =  Wy  ~  <Vy>^  an<^  to  t*ie  sPatial  coordinate  £  (i.e. 
x  we  obtain  the  transformed  equations  of  motion, 

»» 

V  (O 

v  =  Sh(Ou  +  f  <u2>  -  U2)  ;  u  =  -  Qv  ,  (12) 

x  y  22  v  y  y)  y  x  ’  '  ' 

3 

where  we  have  neglected  terms  of  0(s  )  and  higher.  These  equations  lead  to 
(see  Appendix  I), 


v  =  w.  Sin*  -  -=-/0  Sin2*  , 
x  1  3/2 

6ti 


(13a) 


w  . 

u  =  -  Cos*  -  Cos2* 

y 


i2n 


(13b) 


and 


=  -  —  Jcos  fJTT Or  +  $1  -  Costl  +  Icos  f2($./l)SH)]-Cos2$l,  (14a) 

J7al  ^  '  I  12nM  ^  1  I 


X’  -X 


y  -y  =  ^|sin(*+/r|2T)  -  Sin$|  -  — jsin  ^2(  $+/riOT)  j-  Sin2*|  +  <v^>T,(14b) 


Z  -Z  =  V  T  , 
z  ’ 


(14c) 


where  w±2=vx2+h(  S)u2+VE  ( £)  {<uy2>uy-uy3/3} /S2 ,  and  $=VE  (E,)v±2/Q2.  Also 

*=/l"|2t  +  *  where  ♦  is  the  velocity  space  angle  at  t=0.  The  oscillatory  terms 

of  the  order  w  in  the  orbits  are  not  important  except  in  the  derivation  of 

the  Jacobian  of  transformation  from  the  integration  variables  (x,v  ,v  )  to 

x  y 

(E.,w.,$)  which  will  be  necessary  in  the  following.  For  simplicity 
therefore,  we  shall  ignore  the  oscillatory  terms  of  0(w)  everywhere  except 
in  the  derivation  of  the  Jacobian.  This  restriction  can  easily  be  relaxed. 
Using  (13)  and  (14)  in  (8)  we  obtain 


f,(x,v)  =  -e 


;£oa’H)[l' 


dk/exp{ik'xx}<Kk/ ) 


R  v 

r  VJn'(CT,)  (w  ■  kV  )exp(i  [n'(*  -  a'  )  +  k'  l  -  r£-^Sin*|] 

-Jdkx  *<k;>L~ - y*  f - 11 


^  w  -1^ n' 


8-k  <  v  >  -  k  v  1 

y  y  z  z) 


where  a'  =  k.w./Q,  k.2  =  k  2/h  +  k2/t|^,  a  =  tan  3  (k  /t)/k  ) ,  £,  =  x+u  /Q 

1  1  ’  1  x  y  '  x  y  y 

and  J  ,  are  Bessel  functions.  The  projection  of  (15)  in  k  is 
n'  r  j  x 
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fiu  (V)  =  * 


Ik.  v- 
x 


jdkx  *<kx>  Jdx  f0  [  exP{-i<kx  -  kx)xl 


^  J  <,')  Jn(a)(W-  kV  )  ( 

~>  - £expji^(k/  -  k.^)?,  +  (n'  -  n)t+n«-n'a' 

n'n  [ay-JrTn'Q-ky<vy>-kzvz  ) 


The  perturbed  density  is  then  obtained  by  integrating  the  perturbed 
distribution  function  over  all  velocities 


nlk  =  -  ^  jdvxdvydVJdXIdkxf<kx)f0  exp(i(k'-kx)x)-exp(i(k'-kx>5)  F  ,  (17) 


where, 


V  Jn'(a')Jn<ff>  f  ) 

F  =  ( co—k  V  )  )  - — -  exp  [i(n'-n)*+na-n'a' 

y  g  V*  ctt-n'  /t|Q-k  <v  >-k  vK  ' 

n  '  n  1  tt  v  7  7 


y  y  z  z 


Equation  (17)  is  the  most  general  form  for  the  perturbed  density  for  either 
ions  or  the  electrons.  In  the  quasineutral  limit  the  most  general 
dispersion  relation  for  the  electrostatic  waves,  in  the  form  of  an  integral 
equation,  is  given  by  Z  Jexp(ik  x)na1L.dk  =  0,  where  a  represents  the 

(X  X  X  K.  X 


species . 


Now  we  transform  the  integration  variables  in  (17)  from  (x,v^,vy)  to 
($,w^,#)  using  the  appropriate  Jacobian,  which  in  this  case  is  (see 

Appendix  II)  J=-/qw^;  and  expand  the  exponentials  in  x  and  £  around  £.  and 

2 

retain  terms  up  to  0(e  )  to  obtain. 


•'f.rv^ur«#v  *»* 


If 

£ 

:* 


I 

Cv 

iff 

$ 


m 


NeB 
n, ,  =  — 
lk  m 
x 


+  ®  “  2n  +  ®  [_o(  2  2\/9) 

:  J«  \i\  K  JOWS-!  K«ki)  V"2 

»  o  o 


where, 


(19) 


"  2 
Vp(«wf 


G  =  (1-F)  l-i(k'-k  )  ■  v,  —  -i  (k  ' -k  ) 
v  '  V  v  x  x'  4  2^3  )  x  x' 


w^Cos# 

~W~ 


(k'-k  ) 
x  x 


2  2 
2 Cos  $ 

2n»2 


(20) 


We  shall  first  obtain  the  differential  equation  for  the  fluid  K-H  modes  as 
the  fluid  limit  of  (19). 

(a)  Low  Temperature  But  Arbitrary  Shear. 

Consider  k^O  and  en=0.  For  a  quasineutral  plasma  the  condition  that 

the  net  perturbed  density  (19)  vanishes  provides  the  electrostatic 

dispersion  relation.  Further,  since  k(1=0  the  electron  contribution  to  the 

net  perturbed  density  is  ignorable  so  that  the  ion  perturbed  density  set 

equal  to  zero  leads  to  the  desired  general  dispersion  relation.  We  shall 

now  proceed  to  specialize  the  general  dispersion  relation  for  low 

temperature  (k^p^->  0).  Upon  $  integration  the  term  proportional  to  Cos#  in 

G,  vanishes  and  n'  becomes  n  in  (18).  Since  low  temperatures  are  of 

interest  we  keep  n=0,±l  terms  in  F.  Higher  harmonics  are  associated  with 

higher  orders  in  temperature.  The  argument  of  the  Bessel  functions  can  be 

written  as  <j=(w^/vt ) (k^L) t  and,  assuming  that  the  factors  (k^L)  and  (wj/vt) 

are  less  than  or  of  order  unity,  we  can  expand  the  Bessel  functions  to 

9 

0(e  )  so  that, 


8 


V, 

V/tr, 

V>*K 

&£ 


1-F  =  1- 


1  + 


k  ( <v  >-V  ) 
” Z _ g _ 


2  2 
k  wf 

y  -l 
2  2 
2  ns 


1  - 


2  2 
Wj  -  na 


f.,2  ,2  2“l2kxkx  2ikvwlg<kx-kx) 

4ns2  P  x  '  Wl2-  ns2  +  Wl2-  ns2 


(21) 


It  is  important  to  note  that  under  these  conditions  the  terms  of  order 
unity  in  (1-F)  drop  out.  Therefore,  the  term  proportional  to  <Vy>-V^  = 
V  "(w, 2/2S2+p. 2)/2n2  ~  0(e2),  which  is  responsible  for  the  K-H 

Ej  ±  1 

instability,  along  with  the  other  terms  of  the  same  order,  are  the  leading 

terms  of  (19).  This  will  lead  to  the  dispersion  differential  equation 

(22),  describing  the  K-H  modes.  However,  when  a~0(l),  the  Bessel  functions 

2 

can  no  longer  be  expanded  and  consequently  the  role  of  the  0(e)  terms,  and 

in  particular  the  Vg"  term,  will  diminish  in  importance.  This  situation, 

which  corresponds  to  cases  in  which  k^p^>0(l)  (e  may  still  be  small),  may 

lead  to  a  significant  change  in  the  character  of  the  mode,  and  possibly  to 

finite  gyroradius  stabilization  of  the  K-H  instability. 

In  order  to  obtain  the  dispersion  differential  equation  for  the  K-H 

modes  we  need  to  obtain  n^(x)  =  /n^x  exp( ik^xjdk^ ,  for  the  ions  and  set  it 

2 

equal  to  zero.  After  carrying  out  the  w^  integration,  p.  can  be  factored 
out  and  the  resulting  equation  becomes  temperature  independent  in  the  p^-O 
limit.  The  k^  integration  provides  a  delta  function,  2n8(x-£,)>  which  makes 
the  K  integration  trivial  and  converts  the  £.  dependence  into  x  dependence. 
Thus,  performing  the  £,  integration  after  using  Jdk'  x<f>(k^ '  )exp(  ik'  £,)  = 


<p(  £) ,  we  obtain  the  second  order  difletenMa]  equation  for  the  K-H  modes  in 
the  zero  temperature  limit. 


I 


where, 


l  a  2  kvVE(x) 

t  -  *<*>  h  ~  -  ■xi— 


C(x)  )  4>(x)  =  0  , 


2k  VE(x)tt  +  n(x)ST 
A(x)=  «  * 

-  h(x)Q/ 


B(x)=  1  + 


2  2  ’ 

»j-n(x)r 


C(x)  =  -x-i - ■» 

^-n(x>r 


and  «..=w-k  V„(x).  Equation  (22)  is  the  dispersion  differential  equation 
x  y  & 

for  the  fluid  K-H  modes  valid  for  arbitrary  shear  strength  i.e.,  for 
arbitrary  values  of  V  '/S.  as  long  as  h>0.  Starting  from  the  fluid 
equations  one  can  also  derive  (22).  It  should  be  noted  however,  that  when 
n  becomes  small,  the  arguments  of  the  Bessel  functions  become  large.  Thus 
for  h->0  the  expansion  of  the  Bessel  functions  for  large  argument  may  not  be 
a  good  approximation  unless  the  ion  temperatures  and  hence  also  becomes 
vanishingly  small.  Thus  if  h  is  considered  arbitrary  then  the  fluid  theory 
is  valid  only  when  h>0  and  the  arguments  of  the  Bessel  functions  are 
sufficiently  small. 

Pritchett  and  Coroniti^  used  the  fluid  theory  to  obtain  the  dispersion 
equation  for  the  K-H  modes  for  warm  plasmas.  They  incorporated  the 
temperature  effects  through  the  pressure  tensor  and  arrived  at  equation 
(16)  of  reference  (7)  as  the  dispersion  differential  equation.  If  the 
temperature  related  terms  (which  are  included  in  the  factor  R  in  their 
paper)  are  set  equal  to  zero  in  equation  (16)  of.  their  paper^,  then  it 
reduces  to  our  equation  (22).  At  this  stage,  by  defining  a  transformation 
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U=(w1^'-kyS^)/D  where  D=h-(co1/5)z ,  it  can  be  shown  that  all  the  h  related 

terms  can  be  transformed  away  and  (22)  can  be  expressed  in  the  form  of 

(34),  which  represents  the  classical  K-H  mode  except  that  $(x)  is  replaced 

by  U.  Thus,  although  the  eigenfunctions  are  affected,  the  eigenvalues  of 

the  classical  K-H  modes  remain  unaffected  by  the  magnitude  of  the  velocity 

shear  in  the  zero  temperature  limit.  This  is  physically  satisfying,  since 

if  the  gyroradius  is  zero  and  in  the  absence  of  an  equilibrium  x  drift, 

particles  cannot  sample  the  x  direction  and  hence  can  not  experience  the 

velocity  shear  irrespective  of  its  magnitude.  However,  for  finite 

temperatures,  the  magnitude  of  the  velocity  shear  will  play  a  role  as 

evidenced  in  the  simulation  results^  and  the  experimental  results  of 
8 

Jassby  .  Unfortunately  in  the  simple  model  for  the  temperature  assumed  in 
the  reference  (7)  this  cannot  be  explained  since  all  the  temperature 
related  terms  can  also  be  transformed  away'*.  For  arbitrary  shears  the 
temperature  correction  to  the  lowest  order  involves  a  great  many  terms  and 
will  be  discussed  in  a  subsequent  paper.  Here  we  make  the  weak  shear 
approximation  but  treat  temperature  to  be  arbitrary  and  find  that  the 
finite  gyroradius  stabilization  of  the  k-H  modes  can  be  understood  and 
predicted  by  the  kinetic  theory. 


(b)  Weak  Shear  But  Arbitrary  Temperature. 

Now  we  consider  the  limit  where  the  shear  is  not  strong  but  temperature 
and  kn  are  arbitrary.  In  this  limit,  since  the  the  x  dependence  of  the 

t 

equilibrium  quantities  are  weak,  we  can  assume  k  =k  in  (20)  so  that  G=l- 

X  X 

F.  Also  for  weak  shears,  h~l  will  be  considered.  For  generality  we  will 
also  include  an  equilibrium  density  luofNe  N(£)-  The  perturbed  density 
after  carrying  out  the  $  integration  becomes, 


r  rk|  •.!  »«*  < 


4.00  4. 00  4-  00  4. 00 


2ne  0  r  (  \  f  T  ,  ^ 

n(x)=  -  dkxexp[ikxxj  d$  dw^  dvz  dkx' *(kx' >N£S)exp  ji(kx' -kxKj , 

a  J  J 


_co  Q  _  00  _  00 


*  v  ,2, 


.«p(-  f(v^))fT  , 

\  l  l  1  nQ  -k  v  J 

n  1  2a  a  z  z  ' 


Here  a  denotes  the  species  and  «  -  k  V  =  +  «„  -  w  where  w„  = 

y  got  1  2a  a  2a 

M  2  ★ 

k  V_  p  /2  and  00  =  k  e  p  2  .  In  the  denominator  we  have  approximated  w 

y  E  a  ay  na  a  a  vv 

2  2 

-  ky<Vy>  i:  -  a>2,  thereby  replacing  w^  by  its  average  value  2vta  .  Now 

the  w.  and  the  v  integrations  can  easily  be  performed  to  yield, 

_L  Z 


4.  UJ  4.  wu  4-  w 

e  p  f  (  . 

nla(x)=  it*  dkxexP<ikxx>  d^n0(^)  dV*(V)exp[i(kx'~kx)^J 

a  J  J  J 


_  03  —00 


■k 

V  (  «, +w0  -111  1  00.,-w,,  -n2  'l 

y  1  2 * 4“  “  z  y  2°  “  r  (b,  > , 

V  U2|k„|vJ  l  /2  |k||  |vtaJ  “  1 


where  rn(b)=In(b)exp(-b) ,  and  In  are  the  modified  Bessel  functions  and  Z(C) 

2  2  2 

are  the  plasma  dispersion  functions.  Here  b^=pa  (kx  +k^  ).  We  can  expand 

T  in  k  ^p  ^  so  that  T  (b.)  ~  T  (b)+T  (b)k  ,  where  T  =dT  /db 

n  x  a  n  1  n  n  x  a  n  n 

2  2  2 
and  b=k  p  .  Since  we  have  consistently  retained  terms  up  to  e  ,  we 

y 

4  4 

neglect  the  terms  0(kx  )  and  higher  which  are  of  higher  order  in  e. 
Using  this  expansion  for  F^fb^)  in  <?7)  v<=>  perform  the  remaining  intergals 
to  obtain, 


K«W 
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2  2  2 

where  Xa  =mavta  ^^Oofa  ’  *s  tke  De^ye  length  for  the  species  a.  The 
quasineutrality  condition  |e|(n^  -  n^e)=0,  leads  to  the  general  dispersion 

differential  equation  for  the  electrostic  modes  in  the  low  shear  limit. 

2  2 

Retaining  only  the  n=0  term  for  the  electrons  and  considering  p  «p.  we 

G  X 

obtain, 


P2A(X)  ^  +  Q(x)  4>(x)  =0  , 

1  dx^  ) 


where, 


A(x)  = 


y  v 

/2 


Wl+  W2~  W 


lkl|lvi 
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Here  x=l\/Te,  y=nu/me  and  b=(kyp^)  .  The  subscripts  on  and  co  ,  which 
are  for  the  ions,  are  suppressed. 

In  order  to  recover  the  dispersion  differential  equation  for  the  K-H 
modes  in  the  weak  shear  limit  we  set  eo*=0,  k( (=0  and  retain  n=0,+l  in  (30) 
and  (31)  so  that, 


A(x)  = 


Q(x)  =  1  - 


<V  “?]  ' 

— — -  rn(b)  + 

wr  w2>  0 


«,  +  w9 

— - -  I>,(b)  + 

V  “2J  0 


- -  =2 - 2  Fl(b>  ’ 

(Wl-a>2)  -  tf) 


IOj-Wj 


(wr“2)2-  ai 


2rx(b) 


Further,  if  we  take  the  low  temperature  limit  (kyp^  <<1*  thus  ^-(l-b), 

I^-b/2,  T  q  —  1,  T  ^-1/2)  and  low  frequency  limit  /Q^«l)  in  (32)  and 

(33)  and  substitute  these  in  (29),  we  recover  the  starting  differential 

2 

equation  for  the  K-H  modes  widely  used  in  the  literature  , 


,  ,2  9  9  k  V„(x)  9  ) 

2d  ,22  y  E  2  , ,  x 

p.  — ~  -  k  p.  +  — 2 -  p.  >  <f>(x)  =  0  . 
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For  higher  temperatures  (k^p^l)  the  Bessel  functions  can  no  longer  be 

o 

expanded  and  as  explained  earlier  the  terms  of  the  order  e  ,such  as  Uj* 
become  less  important.  Neglecting  Wj  with  respect  to  Mj  and  ignoring  the 
density  gradient  in  (30)  and  (31)  ve  recover  the  starting  equations  for  the 
higher  frequency  ion-cyclotron-like  modes^. 


NUMERICAL  RESULTS 

(i)  Kinetic  Kelvin-Helmholtz  Modes 

We  now  numerically  find  the  eigenvalues  of  the  general  dispersion 
differential  equation  (29),  which  is  still  in  the  weak  shear  limit,  first 
in  the  fluid  limit  to  recover  the  fluid  results  for  the  electric  field 
profile  E(x)=EQtanh(x/L) .  We  then  use  the  parameters  of  the  simulation^ 
to  reproduce  those  elements  of  the  table  (1)  of  Pritchett  and  Coroniti^ 
which  are  accessible  to  the  weak  shear  limit  of  the  theory.  The 
differential  equation  (29)  is  solved  numerically  by  a  shooting  code  for  the 
complex  eigenvalues,  to.  We  assume  WKB  nature  of  the  solutions  for  x-*»  and 
as  the  boundary  condition  demand  that  the  energy  is  outgoing.  With 
and  specified  at  the  boundaries  for  an  initial  eigenvalue  <0q,  we  use 

a  variable  stepsize  integrator  to  obtain  $  and  at  the  origin  where  the 
matching  condition  is  examined.  If  the  matching  condition  is  not  satisfied 
a  new  w  is  assumed  and  the  process  iterates  till  the  eigenvalue  is 
obtained  with  the  desired  accuracy. 


To  recover  the  fluid  limit 
0 


we 


V_=( V_  /v .  )  =  1 ,  u=k.,/k  =0.0001  and  t=3.5. 
E  E  l  '  II  y 

0 


consider  e=0.1,  t  =0,  u=1837, 

n  ’ 

For  kyP^.=0.02,  0.05  and  0.08  we 


obtain  YL/VE  =0. 1369,  0.1868  and  0.1067. 


These  results  coincide  with  the 
fluid  calculations  and  are  consistent  with  the  normalised  growth  rate 
against  the  dimensionless  wavenumhei  pl<>i  for  the  fluid  K-H  instability 
provided  in  figure  (2)  of  reference  (7).  Thus  the  fluid  results  are 
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recovered  from  the  kinetic  dispersion  relation  for  the  parameter  range  that 
are  fluidlike. 

We  now  use  x=l ,  to  match  the  simulation  parameters Also  since  x=l. 

c  =V2v.  where  c  ^=2(T  +T.)/m. .  For  V_=/2V_^/c  =0.764  we  find  that  the 
si  sell  tbs 

eigenvalues  of  (29)  for  kyL=0.393  and  e=0.19  and  0.38  are  given  by  yL/M^= 
0.188  and  0.191  while  the  corresponding  eigenvalues  for  VE=1.513  are  0.184 
and  0.189.  Comparing  these  with  the  corresponding  elements  of  the  table 
(1)  in  Pritchett  and  Coroniti^  we  find  that  our  theory  is  in  agreement  with 
the  simulation  results^.  Figures  (la)  and  (lb)  are  the  eigenfunction  and 
the  profile  for  the  external  electric  field,  for  e=0.19  and  Vg=0.764.  For 
moderate  shears  (e<0.5)  better  agreement  with  the  simulation  can  be 

expected  if  the  assumption  of  weak  shear,  i.e.  kx'=:kx,  is  relaxed.  This 
will  be  the  subject  of  a  future  article.  Higher  shear  values  however,  are 
inaccessible  to  the  theory  at  the  present  simplified  differential  equation 
level.  For  very  high  shears  a  dispersion  relation  in  the  form  of  an 

integral  equation  as  provided  in  (17),  will  become  necessary.  Also, 
uncertainty  in  the  simulation  results  is  expected  for  higher  shears  where  q 
differs  significantly  from  unity  (but  still  positive),  unless  the  initial 
loading  is  in  accordance  with  the  equilibrium  distribution  given  in  (4). 

_  t 

When  shear  is  weak  and  a^V^e  (the  peak  value  of  V^,  (x)/2)  is  small,  q  is 
close  to  unity.  For  low  temperatures  (eO),  the  equilibrium  distribution 
function  (4)  can  be  reduced  to  a  Maxwellian  shifted  in  the  y  velocity  by 

the  magnitude  of  tl.?  E  X  B  drift  which  is  approximately  Vg(x).  Such  a 

distribution  function  was  used  by  Pritchett  and  Coronoti^  for  the  initial 
loading.  While  acceptable  for  small  and  especially  for  low 

temperatures,  this  method  may  lead  to  significant  relaxation  of  the  initial 


equilibrium  for  large  thereby  affecting  the  equilibrium  parameters 

substantially.  Thus,  for  large  o^,  the  interpretation  of  the  final 

simulation  results^  remain  dubious. 

We  now  study  a  different  electric  field  profile,  E(x)=EQsech  (x/L). 
Once  again,  to  check  the  fluid  limit  we  use  e=0.1,  t=l,  p=1837  and 

u=0.0001.  Figure  (2)  is  a  plot  of  the  growth  rates  and  the  real  frequency 
of  the  K-H  modes  as  obtained  from  equation  (29)  (solid  lines)  against  b. 

Q 

The  dashed  lines  are  the  fluid  results  provided  in  Drazin  and  Howard  in 
their  table  (1)  under  sinuous  mode.  Good  agreement  can  be  expected  for 
e<<l.  For  larger  e  however,  the  <*>2  appearing  in  the  denominator  of  (32) 
and  (33)  can  not  be  treated  as  negligible.  For  larger  e,  the  denominator 
can  influence  the  results  by  enhancing  the  growth  rates.  This  important 
temperature  related  effect  cannot  be  predicted  by  the  fluid  theory 
including  temperature,  as  given  by  Pritchett  and  Coroniti^. 

A  two  dimensional  particle  simulation*®  using  an  appropriate  initial 

distribution  function  to  study  the  ion-cyclotron-like  modes  is  currently  in 

progress.  Ultimately  we  shall  compare  the  K-H  modes  with  the  ion- 

cyclotron-like  modes'*  through  numerical  simulation  as  well  as  through 

theory.  Thus  we  will  use  y=100  which  is  being  used  in  our  simulation*®. 

To  study  the  finite  gyroradius  stabilization  of  the  K-H  modes  we  plot  the 

growth  rates  normalised  by  the  ion  cyclotron  frequency  5^  against  b,  for 

various  values  of  u  in  figure  (3).  We  consider  t=3.5,  e=0.43  and  a  mild 

density  gradient  en=-0.05  centered  around  xnQ=1.33p.  such  that  £^=-0.05  for 

x^-pX  x  <'xng+pi  anc*  0  otherwise  and  as  the  d.c.  electric  field  profile  we 

2 

consider  E(x)=EQsech  (x/L)  (see  figure  (4)).  The  growth  rates  of  the  K-H 
modes  are  expected  to  reduce  due  to  the  density  gradient**  but  with  the 
mild  used  here  this  decrease  vi-  n<>t  ignificant.  From  figure  (3)  we 
see  that  for  a  given  e  the  growth  rates  peak  for  a  particular  b.  The  peak 
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of  the  spectrum  is  localized  for  b~0.16  (k  p^~0.4  and  k  L~0.93)  and  is 

y  y 

maximum  for  transverse  propagation  (u=0.0001).  As  the  obliqueness  u  is 
increased,  there  is  a  sharp  decrease  in  the  peak  value  along  with  the 
narrowing  of  the  spectrum.  Beyond  u=0.0075  the  growth  of  the  K-H  mode  is 
substantially  reduced. 

(ii)  Ion-Cyclotron-Like  Modes 

Now  we  consider  the  case  where  k^p^>l.  As  explained  earlier,  in  this 

domain  we  can  no  longer  expand  the  Bessel  functions  and  consequently  the 
2 

0(e)  terms  responsible  for  the  K-H  modes  play  only  a  minor  role.  For  the 

range,  in  which  we  are  now  interested,  where  k^p^~3,  we  can  neglect  these 

terms  for  convenience.  To  explain  the  ion-cyclotron-like  modes  we  first 

resort  to  the  piecewise  continuous  field  profile^  (see  figure  (5)).  This 

is  an  idealization  of  the  actual  field  profile  and  we  use  it  only  to 

demonstrate  the  principles  involved  and  to  obtain  a  good  starting 

eigenvalue  for  numerically  tracking  the  eigenvalues  for  a  smooth  profile. 

For  a  piecewise  continuous  profile  it  is  trivial  to  derive  the  nonlocal 

dispersion  relation.  Setting  and  w*  equal  to  zeros  in  (30)  and  (31),  we 

use  (29)  as  the  differential  equation  for  the  modes  in  question.  In  the 

region  over  which  the  electric  field  is  localized  (we  shall  refer  to  this 

region  as  region  I)  there  is  a  Doppler  shift  in  the  frequency,  i.e., 

w-»w,=w-k  V.,,  while  outside  this  region  where  the  electric  field  is 
1  y  E  ° 

nonexistent  (region  II)  there  is  no  shift  in  the  frequency.  This  is  the 
essential  feature  distinguishing  the  two  regions.  The  matching  condition 
of  the  logarithmic  derivatives  of  the  solutions  of  (29)  at  the  boundary 
x=L/2,  provides  the  nonlocal  dispersion  lelatinn. 


L* J,*  I,'  I.' 


KItan(Kj./2e)  =  iic^  , 


where  =Q(«j)/A(«j)  and  Kjj  is  identical  to  Kj  if  co^  is  replaced  by  w. 

5 

For  details  we  refer  to  our  earlier  papers  .  The  dispersion  relation  was 
solved"*  for  a  wide  range  of  parameters  to  find  growing  modes  distinct  from 
the  K.-H  modes.  We  first  give  a  physical  description  for  the  origin  of 
these  modes. 

The  dispersion  relation  of  the  electrostatic  ion  Bernstein  modes  is"*, 
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where  k.  ~  0  is  assumed.  The  energy  density  of  these  modes  are, 
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Clearly,  these  are  positive  energy  waves.  Introduction  of  an  uniform 

electric  field  in  the  x  direction  initiates  an  E  X  B  drift  in  the  y 

direction  and  consequently  there  is  a  Doppler  shift  in  the  frequency  i.e., 

w-»a>..=a>-k  V_.  The  energy  density  in  the  presence  of  the  Doppler  shift  is, 
jl  y 

U'  ~  wwja(w^),  which  can  be  negative  provided  <*>^<0.  Now  if  we  consider  the 
localized  field  configuration  as  shown  in  figure  (5),  then  it  is  clear  that 
due  to  the  E  X  B  drift  the  energy  density  in  the  region  I  becomes  negative 
while  it  remains  positive  in  regions  II.  A  nonlocal  wavepacket  can  couple 
these  two  regions  so  that  a  flow  of  energy  from  the  region  I  into  the 


region  II  will  enable  the  wave  t 
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Based  on  this  simple  picture  we 


can  predict  some  gross  features  of  the  instability.  As  for  example,  using 
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the  wave-kinetic  description  it  is  possible  to  obtain  the  energy  balance 
condition  for  the  system  from  which  important  scalings  governing  the  growth 
rate  can  be  predicted.  The  growth  of  the  wave  in  region  I  implies  a  loss 


of  energy  from  that  region.  By  conservation  of  energy,  this  must  be  the 
result  of  convection  of  energy  into  region  II  and  any  local  energy 
dissipation  (S  )  or  free  energy  release  (S+)  processes  in  region  I.  The 
rate  of  growth  of  the  total  energy  deficit  in  region  I  is  proportional  to 
the  growth  rate  y,  the  wave  energy  density  Uj  in  region  I,  and  the  volume 
of  region  I,  represented  here  by  the  extent  in  the  x  direction  (L)  of 
region  I  times  a  unit  area  A^  in  the  plane  perpendicular  to  x.  The  rate  of 
convection  through  A^  is  just  V^U^,  where  is  the  group  velocity  in  the 
x  direction  and  is  the  wave  energy  density  in  region  II.  We  can  then 
write  the  energy  balance  condition  as, 


yUILAi  =  (S+-S_-  VGUII)Ai  , 


where  S+  and  S_  represent  the  source  and  the  sink  in  the  region  I.  The 
eigenvalues  obtained  from  (35)  are  expected  to  satisfy  the  energy  condition 
(38).  For  the  situation  presently  under  consideration  we  do  not  have  any 
external  source  of  free  energy  and  since  k(|~0  the  natural  dampings  are 
absent  and  therefore  S  =S  =0.  Now  it  is  clear  from  (38)  that  if  UT  is 
negative  then  y  can  be  positive  and  hence  the  growth  of  the  wave  is 
sustained  by  convection  of  energy  into  the  region  II  from  the  region  I.  On 
the  other  hand  if  is  positive  then  the  convection  of  energy  out  of  the 
region  I  would  lead  to  a  negative  growth  rate  and  therefore  to  damping  of 
the  waves.  For  S  +=0,  an  important  •■’•dim;  -  in  be  predicted  from  (38)  i.e., 
y/V  *  1/L  which  with  proper  normalizations  can  be  written  as  Im(k  p.)  *  z. 
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In  figure  (6)  we  plot  Im(kxp^)  against  e  and  confirm  this  scaling.  In 
figure  (6)  other  parameters  are;  Vg=2.9,  r=l,  u=0.0001,  u=1837  and  the 
growth  rates  have  been  maximised  in  b. 

We  shall  now  study  the  ion-cyclotron-like  modes  for  smooth  profiles. 
For  this  we  use  an  electric  field  profile  given  by, 


E(X)  =  - ^ -  ,  (39) 

A  sinh  (x/a)+l 

2 

where  A=l/sinh  (xQ/a)  and  Xq=L/2.  For  a-O  (39)  represents  a  "Top  Hat" 
profile  which  reduces  to  half  of  it's  peak  value  at  x=Xq.  As  a  increases 
the  profile  becomes  smoother  and  ultimately  when  a=XQ/sinh-^(l)  which  makes 
A=1 ,  the  expression  (39)  reduces  to  E(x)=EQsech  (x/a).  The  shooting  code 
used  for  the  determination  of  the  eigenvalues  can  operate  best  when  the 
initial  guess  for  the  eigenvalue  is  not  too  far  away  from  the  actual  one. 
Thus  it  becomes  necessary  to  use  (39)  so  that  in  the  limit  a-O  we  have 
excellent  guess  values  obtained  analytically  from  the  nonlocal  dispersion 
relation  (35).  Starting  with  the  eigenvalue  for  the  a-O  case  we  slowly 
increase  a  to  track  the  eigenvalues  for  the  profiles  with  the  desired 
smoothness.  For  b~8  we  have  to  retain  n=0,±l,+2,±3,±4,+5,+6  harmonics  and 
the  associated  plasma  dispersion  functions  in  (29)  which  are  evaluated 
numerically;  thus  the  computations  for  each  eigenvalue  is  CPU  time 
intensive. 

In  figure  (7a)-(7d)  we  display  the  transition  of  the  electric  field 
profile  from  nearly  piecewise  continuous  to  smooth  for  four  different 
values  for  a.  Here  s=0.3.  f  >>  f  ut"i  e  (8a)-(8d)  we  display  the 

corresponding  wavepackets.  Other  parameters  are;  b=8,  x=3.5,  u=0.011, 


U=1837,  V„=0.6,  x  «=1.66p.  and  e  =-0.07  if  x  „-0.9p.<  x  <x  ri+0.9p.  and  zero 


otherwise.  The  growth  rates  did  not  vary  much  during  this  transition.  For 

a=0.1,  0.707,  1.41  and  1.89  the  corresponding  growth  rates  normalised  by 

the  ion  cyclotron  frequency  y/2^=0.048,  0.05,  0.037  and  0.031.  There  is 

only  a  40%  reduction  in  the  growth  rate  from  the  sharp  to  smooth  profile 

and  initially  in  going  from  a=0.1  to  a=0.707  there  is  a  slight  increase  in 

the  growth  rate.  This  is  in  contrast  to  the  K-H  instability  where  the 

growth  rates  are  dependent  on  the  second  derivative  of  the  electric  field 

and  are  therefore  very  sensitive  to  the  scale  size  variation. 

In  figure  (9)  we  provide  a  plot  of  the  growth  rate  and  the  real 

frequency  of  the  ion-cyclotron-like  modes  normalised  by  the  ion 

gyrof requency  as  a  function  of  b.  Here  the  profile  in  (39)  is  used  with 

a=1.87  and  the  rest  of  the  parameters  are  identical  to  figure  (8).  We  find 

that  the  instability  is  peaked  around  kyP.,-3  which  for  e=0.3  corresponds  to 

kyL~10  which  is  an  order  of  magnitude  larger  than  the  corresponding  value 

of  the  peak  for  the  K-H  modes.  Further  for  u=0.011  used  here,  the  K-H 

modes  are  expected  to  be  non-existent  and  thus  the  domain  for  dominance  for 

the  two  modes  can  be  quite  distinct.  This  contradicts  the  conclusion  in 

reference  (12)  where  a  simulation  based  on  only  one  set  of  parameters 

obtained  from  the  idealized  field  profile^  was  used  to  conclude  that  the  K- 

H  mode  will  always  dominate  the  ion-cyclotron-like  modes.  Further  the 

12 

initial  loading  in  the  simulation  (assumed  to  be  similar  to  that  of 

reference  (7))  is  improper  since  for  the  parameters  used  was  extremely 

12 

large  (greater  than  unity),  and  consequently  the  simulation  showed  a 

strong  relaxation  of  the  initial  nonequilibrium  velocity  profile. 

In  figure  (10)  we  use  the  parameters  for  our  simulation^!  to  be 

discussed  in  a  separate  article)  i.f-..  u  =  K>0,  x=3.5,  e=0.43,  u=0.038,  e  =- 

0.05  for  x  A-p.<  x  <x  A+p.  and  "  ntheiwi.se  and  >;  ,  =  l.33p.  to  plot  the 
nO  l  nO  l  nO  l  p 

growth  rate  and  the  real  frequency  normalised  by  the  ion  gyrof requency . 


Here  for  completeness  we  also  include  the  ^  term  in  (29)  to  compute  the 

growth  rates  and  use  exactly  the  same  d.c.  electric  field  and  the  density 

profiles  as  was  used  to  produce  figure  (3).  The  inclusion  of  the  ^  term 

does  not  change  the  eigenvalue  by  much.  The  peak  of  the  spectrum  is  around 

b~14.  In  figure  (3)  we  used  the  same  parameters  to  conclude  that  the 

growth  rate  for  the  K-H  modes  are  reduced  significantly  for  u>0.0075  and 

the  peak  of  the  spectrum  is  around  b~0.2.  Once  again  the  domain  of 

dominance  for  the  K-H  and  the  ion-cyclotron-like  modes  are  quite  distinct. 

Finally  in  figure  (11)  we  provide  a  plot  of  the  real  and  imaginary 

parts  of  the  eigenf requency  to  normalised  by  the  ion  gyrof requency  as  a 

function  of  Vg,  the  peak  value  of  the  equilibrium  E  X  B  drift  velocity 

normalised  by  the  ion  thermal  velocity.  Here  b=10,  t=3.5,  u=1837,  u=0.011, 

e=0.3,  x  „=1.66p.,  e  =-0.07  when  x  „-0.9p.<  x  <x  ,.+0.9 p.  and  0  otherwise. 
nU  l  n  nu  i  no  i 

For  the  external  electric  profile  we  use  (39)  with  a=1.87.  We  see  that  the 
real  frequency  is  almost  linearly  proportional  to  Vg  which  is  in  keeping 

4 

with  the  experimental  results  of  Sato  et  al.  . 


DISCUSSION 

Ue  have  provided  a  kinetic  theory  to  study  the  electrostatic  waves  that 

can  be  excited  in  a  collisionless  magnetized  warm  plasma  by  a  transverse 

velocity  shear.  For  k  p.<<l  we  recover  the  fluid  K-H  modes  and  when  k  p. 

y  l  y  i 

is  increased  we  find  that  the  growth  rates  for  the  K-H  modes  are  reduced 

and  for  large  enough  k^p^  the  K-H  instability  is  completely  damped. 

Further,  the  growth  of  the  K-H  modes  is  severely  affected  by  the  the 

parallel  dynamics.  It  seems  that  for  a  collisionless  plasma  the  K-H  modes 

can  grow  only  for  very  small  As  k..p.  becomes  of  order  unity  the 

expansion  of  the  Bessel  functions  u <  n<>  longer  possible.  Consequently  the 
2 

terms  of  0(e  )  responsible  for  the  K-H  modes  diminish  in  importance.  At 


m 


WA 


this  point  the  large  k^p^  ion-cyclotron-like  modes  dominate.  Further, 
larger  k|(  and  density  gradients  inhibit  the  classical  K-H  wave  growth  while 
both  these  effects  favor  the  ion-cyclotron-like  waves. 

An  important  feature  of  the  ion-cyclotron-like  modes  is  the  fact  that 
the  real  frequency  of  these  waves  are  roughly  around  k^Vg^/2  (see  figures 
(9)  and  (11)).  This  is  similar  to  the  K-H  waves  and  therefore  the  two 
instabilities  cannot  be  distinguished  by  the  scaling  of  the  real  frequency 
with  k  V„. 

y  E 

The  linear  dependence  of  the  real  frequency  of  the  ion-cyclotron-like 

modes  on  the  d.c.  electric  field  was  not  explicitly  discussed  in  our 

previous  papers^.  This  could  have  contributed  to  a  misunderstanding  which 
12 

led  Pritchett  to  conclude  that  since  the  modes  in  his  simulation  for 
pVL=0.3  and  for  k^p^  0.47  and  0.94  displayed  the  linear  dependence  of  the 
real  frequency  on  the  equilibrium  flow  velocity  they  could  not  be  the  ion- 
cyclotron-like  modes  that  we  have  discussed.  A  similar  misunderstanding 
was  also  displayed  by  Sato  et  al.^  in  discussing  their  experimental 
results . 

Since  the  initial  electric  field  profile  used  in  reference  (12)  was  not 
in  equilibrium,  the  system  immediately  relaxed  (see  Appendix-I,  condition 
(A22);  here  a. =( p. /a)V„=(2 . 4)3=7 . 2  >>  1)  to  what  is  shown  in  figure  (4)  in 

II  L 

reference  (12)  which  is  much  different  from  the  initial  profile  given  in 
12 

their  equation  (2).  In  fact  the  initial  profile  is  characterized  by  two 

scale  lengths  L  and  a  with  a  peak  value  of  about  3v^  while  the  final 

2 

relaxed  profile  is  more  like  a  Gaussian  or  a  sech  (x/L)  type  characterized 

by  only  one  scale  length  L  and  with  a  peak  value  of  around  2v..  Also,  as 
12 

explained  the  spatial  extent  I,.  "f  tlic-  electric  field  increased  during 
the  course  of  the  simulation.  ('<>n- ■  i  .if  i  *■  ly  estimating  the  broadening  of 
L  only  by  20%,  and  considering  that  the  final  profile  is  approximately 
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similar  to  sech  (x/L)  with  p^/L=0.25,  then  the  modes  at  kyp^=0.47  and  0.94 
would  correspond  to  kyL=1.88  and  3.76  respectively.  These  values  will  be 
larger  if  the  spatial  extent  is  broadened  more  than  the  20%  assumed.  As  we 
have  shown  through  the  analysis  of  the  kinetic  K-H  modes  and  Drazin  and 

9 

Howard  through  the  analysis  of  the  fluid  K-H  modes  for  shear  profile  of 
2 

the  type  sech  (x/L),  the  K-H  modes  are  strongly  damped  for  k^L>l  and  almost 

non  existent  beyond  kyL-2.  Thus,  the  mode  at  k^p^=0.47  can  be  the  tail  end 

of  the  k-H  spectrum  but  the  mode  at  kyp^=0.94  seems  to  be  completely  out  of 

the  theoretically  predicted  k-H  spectrum  and  the  growth  rate  in  the  second 

mode  is  higher.  Hence,  the  conclusion  based  on  the  simulations  of 

reference  (12)  that  the  K-H  mode  will  always  dominate  over  the  ion- 

cyclotron-like  modes  for  a  configuration  with  a  localized  electric  field 

perpendicular  to  an  external  magnetic  field  is  at  best  dubious.  However, 

12 

we  do  agree  with  the  other  conclusion  that  the  idealized  field  profile 
(piecewise  continuous)  used  earlier^  to  demonstrate  the  physical  principles 
involved  is  not  suitable  for  simulation  purposes  and  that  a  strong 
relaxation  from  such  a  profile  to  a  smoother  profile  is  likely.  The  fact 
that  the  piecewise  continuous  field  profile  is  an  idealization  was 
emphasized  in  our  earlier  paper^.  Here  we  have  provided  an  equilibrium 
distribution  (4),  which  if  properly  loaded  should  ensure  a  good  equilibrium 
even  for  moderate  shears.  Since  the  equilibrium  distribution  as  provided 
in  (4),  is  an  implicit  function  of  £,  the  guiding  center  position,  it  is  not 
in  a  convenient  form  for  initial  loading  in  a  particle  simulation.  For 
this  purpose  we  will  express  (4)  in  terms  of  the  real  position  x.  From  the 
definition  (3)  for  £  we  get, 


ms 


Expanding  V„(S)  around  x  we  can  show  iteratively  that, 


v  -  VE(x)  VE(x) 

s-x '  ‘  V” 


Comparing  (40)  and  (41)  we  find  that, 


Vv  ~  VE(x)  2 

Vy  -  V*>  -  -~'n(x)  *  °<E  >  • 


By  definition  Uy=v^-<v^>  and  using  (42)  along  with  the  expression  (11)  for 

2 

<v  >,  we  can  express  u  =(v  -Vg(x))/h(x)  -  0(e).  Also  expanding  t)( K) 

y  j  j 

2 

around  x  it  can  be  shown  that  tl(  ^)  =  0(x)+0(  e) .  Using  these  to  express  w^ 
provided  in  (6),  in  terms  of  x  we  get, 


k  -  vE<*>y 


+  0(e  ) 


and  therefore  the  equilibrium  distribution  expressed  in  terms  of  x  becomes 


fo=nofoifoll  where  f0ir'^/2neXp£_Pvz  /2)  and 


2rt£0l= 


&  exp 


(  er  2 
r  2  [ 


2  (vy  ^  V*» 


l>h. 


/n(x) 


1  +  0(e)  ). 


It  should  be  noted  that  when  expressed  in  terms  of  x  the  h  dependence  in 
w ^  changes  from  multiplying  the  y  component  of  the  velocity,  to  dividing 
it.  The  distribution  given  in  (43)  is  the  zeroth  order  distribution 
function.  The  correction  to  order  e  is  given  in  Appendix-Ill . 

r 

Consider  the  case  where  oi^Vq/LQ,  the  peak  value  of  the  quantity  /Q, 
is  much  smaller  than  unity  (weak  shear).  Now  ti-»l  and  if  0(e)  corrections 
are  to  be  ignored  then  (43)  becomes  a  Maxwellian  shifted  by  the  magnitude 
of  the  E  X  B  velocity  in  the  y  direction.  Such  a  distribution  was  used  by 
Pritchett  and  Coroniti^  and  is  acceptable  for  weak  shears  (a.<<l) 
especially  for  low  temperatures.  To  find  the  correction  due  to  a^,  we 

t 

express  l/h(x)=l-VE  /S  along  with  the  assumption  that  the  temperature  of 
the  system  will  also  be  affected  so  that  p->P+SfJ  such  that  Sp~0(V_  /Q). 

hi 

t 

Using  these  approximations  and  &P=PVE  /2S  we  can  express  (43)  as, 


£ai-  2!  “■>(-  ![  '«  *  <VV*»2  ])[u^Tf-(<VvE<x))2-  \)  ]  <44> 


The  correction  term  proportional  to  was  also  discussed  in  the  reference 

(7)  but  it  was  not  used  for  the  initial  loading  since  it  was  expected  that 

the  system  would  make  the  necessary  adjustments  and  that  these  would  be 

small,  as  long  as  aj<<l.  Thus  as  long  as  is  small  the  use  of  a  shifted 

Maxwellian  appears  to  be  acceptable,  although  (44)  describes  a  better 

initial  distribution.  For  moderate  shears,  however,  the  particle  loading 

should  be  in  accordance  with  (4t).  otherwise  strong  relaxation  from  the 
12 

initial  profile  will  be  inescapable.  Such  strong  relaxation  from  a 
nonequilibrium  starting  condition  is  in-.-niiahly  accompanied  hv  substantial 
free  energy  release,  which  leads  '>•  .i  dynamic  state  quite  different  from 


the  quiet  equilibrium  essential  for  simulation  of  an  instability.  A 


further  improved  initial  distribution  function  with  the  0(e)  corrections 

included,  is  provided  in  Appendix-III . 

It  should  be  remarked  that  in  most  of  the  experiments^  and  space 

3 

observations  there  exists  a  magnetic  field  aligned  current  in  addition  to 

the  transverse  localized  electric  fields.  In  the  case  of  an  oblique  double 

layer  the  magnetic  field  aligned  current  can  originate  due  to  the  d.c. 

electric  field  component  in  the  direction  of  the  magnetic  field  provided 

there  are  some  collisions.  As  for  example,  in  the  experiments  of  Alport  et 
4 

al.  the  double  layer  has  a  component  in  the  direction  of  the  external 

magnetic  field  which  is  larger  than  the  perpendicular  component,  thereby 

providing  a  large  magnetic  field  aligned  current  also.  Further,  in  some 

13 

recent  space  observations  ion-cyclotron-like  oscillations  have  been 
reported  in  conjunction  with  simultaneous  observation  of  a  perpendicular 
component  of  a  d.c.  electric  field  and  a  magnetic  field  aligned  current  for 
situations  where  the  magnitude  of  the  magnetic  field  aligned  current  is 
below  the  threshold  necessary  for  the  excitation  of  the  current  driven  ion 

14  15 

cyclotron  instability  .  A  recent  study  on  the  effect  of  the 
perpendicular  electric  field  on  the  current  driven  ion  cyclotron 

14 

instability  indicates  that  the  perpendicular  component  of  the  electric 
field  can  lower  the  threshold  for  the  current  driven  ion  cyclotron 
instability.  The  necessary  condition  for  the  current  driven  ion  cyclotron 
instability  is  that  the  parallel  phase  velocity  e*>/k(|,  of  the  ion  cyclotron 
waves  resonate  with  the  parallel  electron  drift  V^,  such  that  (o>-k||V(j)<0. 
For  subcritical  V^,  ( <*>— k | ( ) >0  and  therefore  the  Landau  damping  cannot  be 

14 

inverted  .  For  simplicity  again  consider  the  idealized  field  profile  as 

given  in  figure  (5).  The  introduction  of  the  perpendicular  component  of 

the  electric  field  initiates  a  E  D  diiLt  and  consequently  there  is  a 

Doppler  shift  in  the  frequency  «  i.e.,  w-no..  =w-k  V„  in  the  region  I  over 

i  y  b 


which  the  electric  field  is  localized  in  the  perpendicular  direction.  Now 
the  necessary  condition  for  the  onset  of  the  current  driven  ion  cyclotron 


instability  in  the  region  I  becomes  (w^-k||V^)<0,  which  can  be  satisfied 
even  though  remains  positive.  Thus  the  threshold  value  for  the 

magnetic  field  aligned  drift  necessary  for  the  onset  of  the  current 

driven  ion  cyclotron  instability  is  lowered. 

For  the  cases  where  there  is  a  magnetic  field  aligned  drift  in  addition 

to  the  transverse  localized  electric  field,  the  term  S+  in  (38)  is  not  zero 

and  can  roughly  be  estimated  (using  the  local  theory)  to  be  proportional  to 

LUty,  ,  where  the  local  growth  rate  in  the  region  I,  v1=-QT/Qri  ,  evaluated 
11  II  K(a) 

at  w=w  .  Qd  and  QT  are  the  real  and  the  imaginary  parts  of  the  local 
dispersion  relation  identical  to  the  expression  given  in  (31)  with  and 
w*  set  equal  to  zero,  and  is  the  «  derivative  of  Q^.  In  the  ion  rest 

frame  the  field  aligned  drift  provides  an  additional  Doppler  shift  in 

the  electron  term  so  that  Qj  is  proportional  to  (w^-k^V^).  Assuming  that 
the  field  aligned  current  is  also  localized  within  the  region  I  so  that 
^Rw=^I^Wr  we  can  wr*te  t^ie  energy  balance  condition  as, 


**1»I  -  -  <VVd>“rLV  VCUIIAi  ' 


We  have  neglected  the  ion  Landau  and  cyclotron  dampings.  First  consider 
the  case  where  the  electric  field  is  not  strong  enough  to  make  w^<0  and 
therefore  Uj>0  but  is  less  than  co.  Since  co^<«  it  is  possible  to  have 
(<A)^-kyVg)<0  when  (a>-kyVg)>0  and  hence  the  first  term  in  the  right  hand  side 
of  (45)  provides  a  growth  even  for  subcritical  while  the  convection 
leads  to  damping.  Now  if  <0  and  c'>nM-vin/-n  t  ly  U^<0,  the  convection  will 
lead  to  growth  and  the  first  term  in  'In  >  ight  hand  side  will  contribute  to 
damping.  However,  if  k^V^O  (which  can  be  achieved  by  keeping  constant 
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but  changing  the  direction  of  parallel  propagation  or  vice  versa)  growth  in 


the  region  I  may  be  expected  from  both  the  right  hand  terms  for  («j- 
R| |Vd > >0 .  This  can  be  a  likely  scenario  for  most  of  the  experiments  and 
space  observations  involving  the  ion-cyclotron-like  oscillations  for  an 
equilibrium  that  contains  a  d.c.  electric  field  in  addition  to  a  uniform 
magnetic  field.  More  details  will  be  provided  elsewhere. 


CONCLUSIONS 


Using  a  kinetic  approach  we  have  studied  the  generation  mechanisms  for 

the  electrostatic  waves  in  a  magnetized  warm  plasma  including  a  d.c. 

electric  field  perpendicular  to  the  external  magnetic  field.  Two  distinct 

generation  mechanisms  are  discussed  (i)  Kelvin-Helmholtz  mechanism  and  (ii) 

positive  negative  energy  wave  coupling  mechanism.  The  Kelvin-Helmholtz 

mechanism,  first  discussed  about  a  century  ago1,  depends  directly  on  the 

second  derivative  of  the  d.c.  electric  field  while  the  other  mechanism5, 

depends  on  the  inhomogeneity  in  the  energy  density  of  the  waves.  The  K-H 

instability  can  dominate  for  small  k^p^  if  the  propagation  is  nearly 

perpendicular.  For  a  collisionless  plasma  the  K-H  instability  is  strongly 

damped  even  if  k(|  is  a  tiny  fraction  of  ky.  In  the  theory  we  have  shown 

that  the  terms  responsible  for  the  K-H  wave  growth  are  proportional  to 

"  2 

V_  (x)  and  are  of  order  e  .  Only  when  k  p.<<l  the  Bessel  functions  can  be 
fc  y  1 

expanded  for  large  argument  and  the  order  unity  terms  drop  out  thereby 

2 

making  the  order  e  terms  primary  which  then  assures  the  dominance  of  the 

K-H  instability.  When  k^p.  is  increased  and  is  of  the  order  of  or  greater 

than  unity  the  Bessel  functions  can  no  longer  be  expanded  and  consequently 
2 

the  order  e  terms  responsible  fo:  Mn-  1  H  waves  can  not  gain  prominence. 
At  this  stage  inhomogeneous  energ;.  d<  u  •-  i  diiven  modes’  dominate.  Also, 


the  dominance  of  the  K-H  modes  can  be  reduced  even  for  small  k  p.  if  more 

y  i 


oblique  propagation  (larger  k(|)  are  considered.  Here  ve  have  also  shown 
that  the  inhomogeneous  energy  density  driven  modes  can  tolerate  larger  k|(. 
Thus  the  two  modes  are  quite  distinct  and  depending  on  the  parameter  range 
(system  size,  temperature,  density  gradient  etc.)one  or  the  other  can 
dominate. 

It  should  be  pointed  out  that  while  the  interpretation  of  the 

inhomogeneous  energy  density  driven  modes  is  quite  convincing  for  the  "Top 

Hat"  like  profiles  as  evidenced  in  figure  (6),  it  is  not  so  clear  cut  for 

the  smooth  profiles.  As  the  profile  is  made  smoother  additional  physics  is 

introduced  through  various  resonances  that  are  now  possible  since  can 

now  vary  smoothly  over  a  wide  range  of  values  as  opposed  to  one  constant 

value  in  the  region  I  and  a  different  constant  value  in  the  region  II  for 

the  "Top  Hat"  profile.  Geometry  related  effects  can  also  play  a  role.  It 

was  also  noted  that  as  the  smoothness  of  the  profile  was  increased  it  was 

necessary  to  maintain  a  very  small  amount  of  the  density  gradient  in  the 

transition  zone  in  which  the  electric  field  is  reducing  to  zero,  to 

preserve  the  growth  rates.  This  however,  makes  the  model  more  physical 

4 

since  in  actual  experiments  (e.g.  see  Alport  et  al.  )  a  density  gradient  is 
present  in  the  transition  zone.  It  appears  that  the  density  gradient  acts 
as  a  catalyst  by  maintaining  the  growth  rate  without  much  affecting  the 
real  frequency,  although  the  exact  role  that  the  density  gradient  plays  is 
yet  to  be  fully  appreciated.  The  important  conclusion  however,  is  the  fact 
that  besides  the  K-H  instability  there  is  another  branch  that  can  also  grow 
with  shorter  wavelengths  and  higher  frequency  in  a  plasma  immersed  in  a 
uniform  magnetic  field  with  a  nonuni foim  transverse  electric  field. 

Finally  we  would  like  to  point  <>nt  that  in  the  small  kp.  limit  the 

y  i 

integral  equation  can  be  exactly  t  •  dm  <■  d  •  <>  'he  second  "idri  differential 
equation  (22).  Thus  the  second  order  differential  equation  level  of 
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description  to  study  the  nonlocal  wave  dispersion  properties  employed  in 
this  paper  is  more  accurate  for  the  K-H  modes  than  the  ion-cyclotron-like 
modes  that  grow  for  large  k^p^  For  greater  accuracy  the  eigenvalues  of 
the  integral  equation  which  will  result  by  using  (17)  as  the  perturbed 
density,  must  be  obtained.  This  will  be  the  topic  of  a  future  article. 
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APPENDIX  I 

In  this  appendix  we  provide  the  derivation  of  the  particle  orbits  to 

2 

0(e  ).  The  x  and  y  component  of  the  equation  of  motion  (1),  can  be  written 


*x  =  %  -  WE(X)  ’ 


Vy  =  -Gvx  ,  (A2) 

where  Vg(x)  =  -  cE(x)/Bp.  Expressing  (Al)  in  the  guiding  centre  frame  f, , 

2 

and  retaining  terms  upto  0 ( e  )  we  get, 


,  1  (vv  -  V«)2  "  3  1 

■  «  M)  (v  -  vE(0  -  ^  Ve(E.)  *  0(E3)  . 


We  now  transform  (A2)  and  (A3)  to  a  frame  moving  with  a  velocity  <vy>  in 
the  y  direction,  (i.e.  v^  -»  u^+<v^>)  so  that, 


r 

vx  =  q  n(Ou  *  n(«(<v  >  -  vE(o) 


+  <v  >  -  v. 


.a))2  . 

! _ if _ JF 


VE(0  •  (A4) 


An  expression  for  <v y >  was  given  in  equation  (11)  in  the  text.  Replacing 

v  by  u  +<v  >  in  the  right  hand  side  of  (11),  we  find  that  <v  >  -  V„(£,)  = 

y  y  y  y  E 

"  2  2  3 

V„  (£,)<u  >/2Q  q  +  0(e  ).  Substituting  this  in  (A4)  and  transforming  v  to 

t-  y  y 

u  +<v  >  in  (A2),  we  obtain  the  equations  of  motion  in  the  transformed 

y  y 

2 

frames  to  0( e  )  , 


V  v  . .  I 


'm  '  .-•  ^  ~T*  m  W.  v  n  m  •  m  -  «  -  M  ^  J*  -*  rjf  r^rATT.  W-  w\  Vm  rm  WJ  ^M^Tt^**"** 


Thus , 


ve<Oa‘ 

u  =  Acosf  -  - 9  cos2$  , 

y  i2na)sr 


u  _  ve(^>a 

-p^-  =  /riAsirvf1  -  - sin2$ 

S  6/hST 


(All) 


(A12 ) 


The  constant  A  is  still  undetermined. 


After  multiplying  Equation  (A7)  by  u  it  can  be  written  as, 


/  -2  2  .  3 

,  f  u  0  u  V„(  £,)  .  u 

dll  “2  +  n(C)C  2^~  "  ~2  (  ~3 


i-  -  <u  >u  =  0 

3  y  y  ) 


))■ 


(A13) 


Using  (A6)  we  can  eliminate  from  (A13)  and  define  a  constant  w^  as, 


2  2 

w.  =  v 

1  X 


2  V*>  (  \  ?  \ 

*  na)uy  -  —  (  -  <uy>uy  )  • 


(A14) 


Substituting  uy  and  from  (All)  and  (A12)  into  (A14)  and  retaining  terms 
2 

up,  to  0(t  )  it  can  be  shown  that,  A=w^//f|.  Thus, 


V°wi 

v  =  w.sin$  -  - ■-  sin2$  , 

X  6n(0  ST 


W1  VE(°W1 

U  =  -  C OS$  — : - - — -  cos2i> 

y  /nco  i?n(  q' 


(A15) 


(A16) 
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With  the  velocities  known  it  is  a  simple  matter  to  obtain  the 
positions , 


x '  -  x  = 


v  dx 


cosf'  -  cos$ 


V^VI 


i2na)V 


cos2$' 


cos2$ 


(A18) 


/n(x)a 

Rewriting  =  /r|Ot'+$  =  AiSt+$,  where  T  =  t'-t  and  $  =  /t)Qt  +  $,  we  get, 


x'  -  x 


/n(Oa 


cos 


(Ai(  O 


fix  + 


- — Y-l  cos  2(*+/n(^)8T)  -  cos2$  }  ,  (A19) 

I2n<0  al  v  '  J 

»  2  2 

where  w  =  (^)w^  /Q  .  Similarly  y'-y  can  also  be  obtained.  It  should  be 

noted  that  when  h<0,  the  orbits  become  unstable. 

For  computer  simulations  where  a  distribution  of  particles  is  in 

consideration  we  can  get  an  order  of  magnitude  restriction  necessary  for 

the  stability  of  the  orbit  of  a  typical  particle  (characterized  by  a 

velocity  v^,  the  thermal  velocity  and  a  displacement  p,  the  gyroradius). 
From  (A3)  it  is  clear  that  as  long  as  the  first  term  in  the  right  hand  side 
which  is  of  order  v  ,  is  dominant  the  form  of  (A3)  is, 


(A20) 


X  =  -  2  n(OX  +  corrections 


where  X=(£-x)  =  (v  -V,,(0)/2*  Therefore  the  restoring  force  is  proportional 
y  t, 

to  the  displacement.  This  ensures  periodic  orbits  which  are  stable.  On 

the  other  hand  if  the  second  term  in  the  right  hand  side  of  (A3),  which  is 
2  "  2 

proportional  to  v  V£  (£)/2  ,  dominates  then  (A3)  is  of  the  form, 


vEa> 


x2  + 


corrections 


(A21) 


The  restoring  force  is  now  proportional  to  the  square  of  the  displacement. 

Hence,  the  orbits  are  no  longer  periodic  and  therefore  become  unstable. 

Thus  as  long  as  the  second  term  of  (A3)  remain  smaller  than  the  first  term 
2  "  2 

i.e.,  v„  V„  (0/20  <h(Ov„,  we  can  expect  stable  orbits  for  the  typical 

t  fci  t 

particles.  This  restriction  leads  to  the  condition  for  stable  orbits  in  a 
simulation, 


f2vtn(0' 


1/2 


VO 


where  R(0  is  the  local  radius  of  curvature 

Simplifying  (A22)  by  replacing  the  guiding  center 

position  x  and  considering  the  electric  field 

Vc,(x)=V  ^f(x/L),  we  can  define  H(x)  such  that, 
h  b 


(A22) 

(=(|vE(0/vE"(0l)1/2). 

position  £,  by  the  real 
profiles  of  the  form 


*•  „*» 


where  x=x/L  and  e=p/L.  In  order  to  have  stable  orbits  so  as  to  avoid  (or 
minimize)  relaxation  of  the  initial  electric  field  profile  used  in  a 
computer  simulation,  H(x)  should  be  positive  for  all  x.  If  this  condition 
is  violated  then  the  orbits  will  become  unstable  and  the  profile  will  relax 
until  (A22)  is  satisfied. 


APPENDIX  II 


In  this  appendix  we  evaluate  the  Jacobian  of  transformation  from  the 

coordinates  (x,v  ,v  )  to  the  coordinates,  (£,w.,*).  This  Jacobian  can  be 
x  y  -L 

written  as, 


J 


3x 

3£ 


3£ 


3v 

at 


3x 

3x 

3w^ 

l,* 

3$ 

3v 

3v 

_ X 

X 

Wj_,  * 

3w^ 

K,* 

3$ 

3v 

3v 

_Z 

_Z 

Wj_,  $ 

3Wj^ 

K,* 

3$ 

,1 


(Bl) 


Using  the  definition  of  £  as  given  in  equation  (3)  in  the  text  we  can 
evaluate  the  elements  of  the  first  row  so  that, 


i  3v 

i  3v 

_  I  _Z 

i 

.  1  — SL 

1  _Z 

Q  3£ 

R  3w^ 

R  3$ 

3v 

3v 

3v 

X 

_ X 

X 

3£ 

3w^ 

3$ 

3v 

3v 

3v 

__Z 

_ Y. 

_Z 

3£ 

3w^ 

3$ 

(B2) 


where  we  have  suppressed  the  subscripts.  Multiplying  the  last  row  by  1/Q 
and  adding  it  to  the  first  row  we  get. 


Thus  the  determinant  has  been  considerably  simplified  and  can  be  expanded 


J  =  HU) 


3v  3v  3v  3v 
x  y  x 


3w^  3*  3$  3w^ 


Recall  that  v  =u  +<v  >  and  using  the  expressions  for  v  and  u  from 

y  y  y  x  y 

(A15)  and  (A16),  and  using  equation  (11)  from  the  text  for  <vy>  we  get, 


II  o 

VEa)wl 

v  =  w.sin*  -  - 7  sin2*  , 

6n(S)J/'V 


W1  (  cos2* 

/n(0  E  4n(^)2az  ^  3 


)  • 


The  derivatives  necessary  in  (B4)  can  be  easily  obtained,  and  retaining 
2 

terms  up  to  0(e  )  we  get, 


'x  VEa)wJ 

—  =  sin*  -  — - - 

*i  3n< 


S 


In  this  appendix  we  provide  an  equilibrium  distribution  function 

suitable  for  initial  loading  for  a  particle  simulation  studing  the 

electrostatic  waves  due  to  an  equilibrium  field  configuration  containing  an 

uniform  magnetic  field  and  a  perpendicular  component  of  a  nonuniform 

electric  field.  We  shall  include  the  corrections  which  are  of  the  order  e 

2 

but  ignore  the  order  e  corrections  for  the  time  being.  Expanding  q( £) 
around  x  we  get  q( £)=h(x)+( S-x)n' (x)+0( e2) ,  where  ty=VE  /».  Using  (41) 
for  (£-x)  we  can  write  down, 


n<5)  =  n<x)  +  ( 


Vy  -  VE(X)  )  VE(X) 


n(x)a 


)^  +  o(e2)  . 

'  r\ 


Als4  from  the  expression  for  uy  given  in  (A16)  we  find  that  the  time 

2  2  2  2 
average  <uy  >=w^^  /2q(S)-  Expressed  in  terms  of  x,  <uy  >=w1  (l-(Vy- 

Vg(x) )V  E(x)/r|2(x)S22)/2n(x)+0(e2) .  Using  this  and  (Cl)  in  the  expression 
2 

for  w^  as  provided  in  (6)  we  get, 


.2  -  v2  .  b  -  V,))\  L  .  v  (x)H.i  .  (vy  -  VE(X>)21  V!! 

W1  vx  +  n(x)  l  y  Vx,Jl2  2h(x)  J  „2 ,  . 


2n(x)  J  2  2  + 

m  '  1  ry (x)sr 


Using  (C2)  in  the  expression  for  f-.  and  expanding  the  0(e)  terms  we  get, 


gf  2  <yv*»4 

2  vx  +  n(x) 


/q(x) 


BMW 


m 


Fig.  2  —  The  real  and  imaginary  frequencies  normalized  by  ft,  for  the  K -H  instabilities  for  the  d.c. 
electric  field  profile  given  by  £  =  £0  sec/t 2  (x  /L )  are  plotted  as  a  function  of  b .  The  solid  lines  are 
the  eigenvalues  of  the  equation  (29)  while  the  dotted  lines  are  the  fluid  results  given  in  the  reference 
(9).  Here  t  =  0. 1,  r  =  1,  n  =  1837  and  u  =  0.0001. 
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Fig.  8  (Continued)  —  The  real  and  imaginary  parts  of  the 
corresponding  eigenfunctions  for  the  profiles  in 
figure  (7).  Here  b  =  8.  ^  =  1837,  Vt  =  0.6.  t 
=  3.5,  u  =  0.01 1  and  (n  —  —0.07  if 
-*„(»  -  0.9p,  <  r  <  .v,l0  +  0.9P,  and  x„„  = 

1.66„,  .(c )a  =  1.41  and  (d)  a  =  1.89. 


Fig.  1 1  —  The  normalized  real  and  imaginary  parts  of  the  eigenfrequency  of  the  ion-cyclotron-like 
waves  plotted  against  VE.  Here  b  =  10,  r  =  3.5,  n  =  1837,  u  =  0.011,  e  =  0.3,  a  =  1.87  and 
f„  =  -0.07  if  xn0  -  0.9 p,  <  x  <  xn0  +  0.9 p,  and  zero  otherwise. 
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